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Abstract. We propose a renormalization scheme that can be used as a reliable 
method to identify universality classes of absorbing phase transitions. Following the 
spirit of Wilson's block-spin renormalization group, the lattice is divided into blocks, 
assigning to them an effective state by a suitable Boolean function of the interior 
degrees of freedom. The effective states of adjacent blocks form certain patterns which 
are shown to occur with universal probability ratios if the underlying process is critical. 
Measuring these probability ratios in the limit of large block sizes one obtains a set of 
universal numbers as an individual fingerprint for each universality class. 



1. Introduction 

Scale invariance is known to be an important cornerstone for our understanding of 
critical phenomena and provides the basis of renormalization group theory. On a lattice 
a scale transformation can be carried out by coarse-graining the elementary degrees of 
freedom. One of the best known examples is Wilson's block-spin renormalization of 
the Ising model, where several spins are grouped into blocks P-[3]. If the Hamiltonian 
in terms of these block variables has the same form as the original one (possibly with 
different parameters), this procedure maps the system onto itself, forming a group of 
renormalization transformations. Critical phenomena are typically associated with fixed 
points under renormalization group transformations. 

Non-equilibrium systems differ from their equilibrium counterparts in so far as they 
involve time as an additional dimension on equal footing with the spatial degrees of 
freedom. Nevertheless the situation is similar, i.e. by appropriately coarse-graining the 
elementary degrees of freedom we can identify critical phenomena as fixed points under 
renormalization group transformations. In the past decades various renormalization 
techniques for nonequilibrium systems have been developed, including field-theoretic 



renormalization methods density matrix renormalization [5], and continuous non- 
perturbative renormalization techniques [B], to name only a few. 

In the present paper we consider a discrete renormalization scheme in the spirit of 
Wilson's block-spin method that is suitable for systems with continuous nonequilibrium 
phase transitions into absorbing states [THTU]. The most important example of such 
systems is directed percolation (DP) [11], a universality class of absorbing phase 
transitions which plays a paradigmatic role similar to the Ising class in equilibrium 
statistical mechanics. Other examples include the class of compact directed percolation, 
the parity-conserving universality class, voter-type phase transitions, the pair-contact 
process with diffusion, and the Manna class (for a review see e.g. [ID])- Universal 
behavior is also observed at tricritical points [121 [13] and in models with algebraic 
long-range interactions [2]. Recently absorbing phase transitions attracted increasing 
attention since the critical behavior of directed percolation could be reproduced 
experimentally for the first time in turbulent liquid crystals |15] . 

Before going into detail let us briefly sketch how the proposed method for 
renormalization works. Restricting for simplicity to the 1+1-dimensional case, it involves 
the following steps: 

- Consider a sufficiently large system on a chain with L sites and periodic boundary 
conditions starting with a homogeneously active state at criticality. 

- Read out the configuration of active sites at time t and divide the chain into blocks 
of equal size h. A block containing at least one active site is called active (1) and 
otherwise inactive (0), giving a string of L/h bits. 

- Scan the bit string cyclically in order to estimate the probability Pc{b,t) to find 
n randomly selected adjacent blocks with the bit pattern c. For example, three 
adjacent blocks can have eight different patterns c G {000, 001, 010, . . . , 111}. 

This procedure is carried out at different times t and for different block sizes b, averaging 
over many runs to reduce statistical fluctuations. Note that the underlying process is 
not affected, all what is needed is to read out the configurations in regular time intervals 
and to process this data in a separate procedure in order to estimate Pc{b,t). 

The probabilities Pc{b, t) will depend on both the block size b and the actual time t. 
However, as will be shown below, the ratios of these probabilities saturate and can be 
extrapolated to t — )■ oo. These block-size-dependent extrapolated probability ratios in 
turn converge to well-defined universal values in the limit b ^ oo, providing a whole set 
of universal numbers which can be viewed as a fingerprint of the underlying universality 
class. With these numbers denoted as Sc, which are as robust as critical exponents, it 
is possible to identify a universality class in a reliable way. 
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2. Block renormalization scheme for systems with absorbing states 
Symmetry as the guiding principle 

The purpose of any block renormalization scheme is to organize the elementary degrees 
of freedom (the lattice sites) into equal groups and to replace the individual states inside 
each group by a single effective state of the same kind. This coarse-graining has to be 
done in such a way that the symmetries of the model are preserved. For example, in 
Wilson's block renormalization scheme for the one-dimensional Ising model the spins 
are grouped in pairs. One of the two spins is taken as the effective spin of the block 
while the other spin is integrated out and absorbed into the effective interaction of two 
neighboring blocks. Obviously, this procedure preserves the global Z2-symmetry of the 
Ising model. 

Devising a similar method for models with absorbing state transitions, the coarse- 
graining procedure should respect the most salient feature of such models, namely, the 
existence of an absorbing state. An absorbing state is a configuration that can be 
reached but not be left and therefore plays the role of a dynamical trap. Consequently 
any renormalization procedure should map a locally absorbing configuration of lattice 
sites onto an absorbing effective state of the block and conversely every non-absorbing 
configuration onto an active effective state. In other words, the renormalization proce- 
dure has to preserve the property of 'being absorbing'. 

Two-state models: Renormalization by a logical OR 

Many models with absorbing phase transitions are two-state models, where each lattice 
site i carries a binary bit Si G {0,1}, denoting inactive and active sites, respectively. 
The renormalization procedure outlined above groups these bits into equal-sized blocks, 
assigning to each block an effective state in form of a single bit which is again either 
or 1. Therefore, the desired map has to be a Boolean function. 

Since the renormalization procedure has to preserve the absorbing state, a locally 
absorbing configuration 00 ... has to be mapped to 0. Conversely, any non-absorbing 
configuration, where the block contains at least one active site, has to be mapped to 
1. Obviously, the only Boolean function with this property is a logical OR. Thus we 
can conclude that any block renormalization scheme for two-state systems with a single 
absorbing state should be based on a logical OR of the individual degrees of freedom. 

The coarse-graining procedure by a logical OR can be implemented in different 
ways. In the following we will use it to renormalize spatial configuration (snapshots) of 
the system at a certain instance of time, as sketched in Fig. [T] After renormalization 
one obtains coarse-grained bit patterns of adjacent block variables. The aim is to ana- 
lyze their probability distribution, to study possible universal properties, and to use the 
results in order to identify the underlying universality class of the process. 
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Figure 1. Block renormalization of two-state models with absorbing states. The 
figure shows a contact process at criticality starting with a fully occupied lattice. The 
upper part shows the temporal evolution, where each row of pixels represents the 
configuration after 10 Monte Carlo sweeps. At a certain point we read out the spatial 
configuration of actives sites which can be considered as a string of bits. The figure 
shows two subsequent renormalization steps symbolized by blue triangles, each dividing 
the chain into groups of six bits and mapping them onto a single bit by a logical OR. 
Obviously these two steps could be replaced by a single one, grouping 36 bits together, 
demonstrating that these operations form a commutative and associative group. 



Scaling properties under coarse graining 

The critical properties of absorbing phase transitions are well described by a 
phenomenological scaling theory [TU] . Starting point is the postulate that such systems 
are invariant under anisotropic scale transformations 

x^Kx, t^AH, p^A-f^/'^p, A^A-'/^^A, (1) 

where A > is a spatial scale factor, p denotes the coarse-grained density of active sites, 
and A = A — Ac denotes the distance from criticality in the control parameter A. Apart 
from few exceptions [HlllG], time and space scale differently, expressed by a dynamical 
exponent z = i'\\/i'± > 1. Thus the scaling behavior of such a system is determined by 
three independent critical exponents B, u^, v\\ from which all other bulk exponents can 
be derived by simple scaling relation4i|. 

X In principle there could be a fourth independent exponent /?' which is associated with the conjugate 
order parameter. However, in most cases a time reversal symmetry implies /3 = /?', see e.g. Ref. [TU] . 
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By grouping b lattice sites into blocks the lattice spacing a increases by a — )■ ab. 
This means that the spatial coordinate x in Eq. ([T]), which measures the distance in 
units of the lattice spacing, decreases by x x/b. Therefore, the scale factor of each 
renormalization step is A = 1/6 so that the quantities discussed above are expected to 
scale as 

x^x/b, t^b-H, p-^b^/^^p, A^b^Z-'^A. (2) 

Note that the density of active blocks increases under coarse-graining, reflecting the fact 
that a logical OR can only increase the density of active blocks, as illustrated in Fig. [T] 
Moreover, the last equation A — )■ A tells us that the distance from criticality effec- 
tively increases under renormalization, driving the system away from criticality. 



Scaling properties of the probability distribution Pc{b,t) 

Suppose that we have read out the configuration of the system at time t and arranged 
it into blocks of size b as shown in Fig. [1] Scanning over the resulting bit pattern 
(and averaging over several runs) we can now estimate the probability Pc{b,t) to find n 
randomly selected adjacent blocks in the configuration c = {ciC2 . . . c„}. By definition, 
these probabilities are normalized by Pc{b,t) = 1. 

How do these probabilities decay asymptotically as functions of time? To answer 
this question we first note that Pi(l,t) is just the density of active sites on the original 
lattice, hence this quantity is known to decay as 

Pi(l,t) = p(t)~t-°, (3) 

where a = P/u^. If the renormalization scheme outlined above works properly, the 
same should hold for all single-bit quantities after coarse-graining, i.e. we expect them 
to decay asymptotically as 

P,{b,t)r^t-" (4) 

for arbitrary block sizes b, differing only by a 6-dependent proportionality factor in front 
of the power law. This observation can be confirmed easily by numerical simulations. 

Let us now turn to the two-bit probabilities. On the one hand. Poo, Pqi, Pio, and 
Pii are restricted by the normalization condition and the left-right symmetry 

Poo{b, t) + Pot{b, t) + P,o{b, t) + Pii(6, t) = 1 , (5) 

Poi{b,t) = P,o{b,t) . (6) 

On the other hand, the renormalization mapping by a logical OR can be expressed as 

Poi{h, t) + P,o{b, t) + Pn{b, t) = A (26, t) ~ t"" . (7) 

Finally, tracing out the second bit one obtains the relation 

Pio(6,t) + Pii(6,t) = Pi{b,t) ~ (8) 
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Figure 2. Temporal decay of the four-bit probabilities Pc{b,t) for two different 
block sizes b ~ 2 and 6 = 8 in the contact process. Note that all quantities decay 
asymptotically with the same power law. 



with the same power law as in ([7]) but with a smaller proportionality factor. These 
four equations imply that the two-bit probabilities with non-vanishing bit patterns (i.e. 
Poi{b,t), PiQ{b,t), and Pii{b,t)) decay asymptotically in the same way as the density of 
active sites p{t), but with different prefactors in front of the power law. 

Similarly all n-bit probabilities for arbitrary block sizes b are found to decay 
asymptotically in the same way as the density of active sites, i.e., 

Pc(^t)~t"" if c 7^ 000... 0. (9) 

This is demonstrated in Fig. |2] in the case of DP, where we plotted the decay of the 
four-bit probabilities for two different block sizes. Note that this common type of al- 
gebraic decay sets in only when the correlation length of the process exceeds the size 
of the block sequence, i.e. ^± ^ nb. This requires a simulation time t ^ {nb)^, where 
z = I'w/i^i. is the dynamical exponent of the process under consideration. 



Time-independent quotients in the limit t — > oo 

Since all probabilities Pc{b,t) with c 7^ 00 ... decay asymptotically in the same way, 
any quotient of these quantities will saturate at some constant value in the limit t — >■ 00. 
As a convenient choice we decided to study the ratios 

^c(M) := ^^^^^ = (c^OO...O), (10) 

where we used the normalization condition in the denominator. These quantities can 
be interpreted as the conditional probability to find n subsequent blocks of size b in the 
bit pattern c, given that at least one of the blocks is active. By definition, they are 
normalized by 

E S.ib,t) = l. (11) 

c 7^ 00.. .0 
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Figure 3. Main steps of tlie metliod proposed in tiie present worlc. 



As discussed before, the conditional probabilities tend to a constant 



S,{b) := \imS,{b,t) 



(12) 



in the limit t — oo. In this way we have defined stationary block probabilities in a 
non-stationary process. 

Universality of the conditional probabilities in the limit b ^ oo 

Since the conditional probabilities Sc{b) can be viewed as some kind of short-range 
correlation functions, we expect them to be non-universal, refiecting the microscopic 
details of the dynamics. However, the main message of this paper is the observation 
that the block renormalization procedure introduced above drives the conditional 
probabilities towards certain universal values which depend only on the large-scale 
critical behavior of the process under consideration. This means that the conditional 
probabilities converge to well-defined universal values 



in the limit of large block sizes. Note that this limit has to be carried out after taking 
t — )■ oo since the two limits fll2p and fll3|] do not commute. 

As will be demonstrated in the following section, it is possible to carry out the two- 
fold non-commuting limit in numerical simulations which allows one to estimate the 
numbers Sc, following the sequence of steps shown in Fig. [31 As mentioned before, these 
numbers can be viewed as a fingerprint of the underlying universality class. Since our 
method can be applied to any models with absorbing states, we can use these numbers 
to verify whether a given process belongs to a certain universality class or not. 

Exact scaling relations for the two-bit probabilities 

As will be shown now, the universal two-bit probabilities Soi, S^o, are special in so 
far as they can be expressed exactly in terms of the critical exponents. To see this 
recall that -Pi(l) is the density of active sites, while Pi{b) is the renormalized density of 
active sites, defined as the probability to find at least one active site in a block of size 
b. According to Eq. this quantity is expected to scale as 



Sc 



lim Sc{b) 



(13) 




(14) 
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as 6 — 7- oo. This asymptotic power law implies the identity 

lini = 2~l"^- . (15) 

Using the exact relations Pi (26) = 1 - Po{2b) = 1 - Poo(&) and Pi (6) = Poi(6) + Pii(6) 
this identity can be rewritten as 

hm = So, + Sn = 2-^/- . (16) 

b^oo l-Poo(6) 

Together with the normalization condition and spatial reflection symmetry we obtain a 
system of three equations 

Soi + Sn = 2-"/"^ (17) 

Sqi + SiQ + 5*11 = 1 (18) 

Soi = 5io (19) 
with the exact solution 

Soi = S,o = l- 2-^/''^ (20) 
5ii = -1 + 2^-'^/'^^. (21) 

If the critical exponents are known, this allows us to predict the two-bit quantities 
and to verify the accuracy of the method in a numerical simulation. In the remaining 
part of this article we use Monte Carlo simulations to do the same for two well known 
universality classes of absorbing phase transitions, namely directed percolation and 
parity conserving class. 

Note that the conditional probabilities with three or more bits cannot be calculated 
in this way since the scaling laws do not provide enough equations. This is good news 
since these probabilities provide additional information in form of universal numbers 
which are independent of the usual critical exponents. From the numerical simulations 
we also estimate these universal numbers for the above mentioned classes. 



3. Directed percolation 

Verification of the two-bit probabilities 

For directed percolation in one dimension the best known estimates of the critical 
exponents are [17] 

(3 = 0.276486(8) , = 1.096854(4) , i^,, = 1.7733847(6). (22) 



Inserting these estimates into Eq. f l20|) and fl2T|) we can compute the conditional two-bit 
probabilities 

5oi = 5io ~ 0.160310(5) , 5ii ~ 0.67938(1) . (23) 
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Figure 4. Two-bit probabilities for directed percolation. Left: Estimation of the 
conditional two-bit probabilities Soi{b) and Snib) in the limit t ^ oo according to (jl2p . 
Right: Extrapolation of the universal values Soi and 5*1 1 for 6 — > cxd according to (jl3p . 
As can be seen, the values converge to the predicted values in Eq. (^5]) which are 
indicated by the arrows. 



Let us now demonstrate how these values can be obtained numerically by performing the 
two non-commuting limits t -> oo and b ^ oo. To this end we simulate an ordinary bond 
DP process starting with a fully occupied lattice over 10^ Monte Carlo sweeps. In certain 
intervals we read out the actual configuration and apply the suggested renormalization 
scheme, reducing b bits to a single bit by a logical OR operation. Scanning over 
the resulting bit string and averaging over many runs, we estimate the conditional 
probabilities Sc{b, t) as a function of time for the block sizes b = 1, 2, 4, 8, 16, and 32. 

As shown in the left panel of Fig. |H these probabilities tend to certain asymptotic 
values Scip) as t is increased. It turns out that the difference Sdb.t) — Sdb) decays 
asymptotically as - a numerical observation for which we do not yet have an 
explanation. Hence, by plotting Sc{b,t) versus 1/t we can easily extrapolate the curves 
to the limit t — )■ oo. 

The extrapolated values for Sc{b) are then plotted in the right panel of Fig. |H In 
order to perform the second limit 6 — ?■ oo it is again useful to plot these values versus 
1/b and to estimate the limit Sc = \imb^ooSc{b) by linear extrapolation. As can be 
seen, this procedure accurately reproduces the predicted values in Eq. fl2^ which are 
indicated by arrows in the figure. 

Numerical confirmation of universality 

In order to demonstrate the universal character of the extrapolated conditional 
probabilities Sc in the case of bit patterns with more than two bits, we measured the 
3-bit quantities for various frequently used models in the directed percolation class, 
namely the contact process with random sequential updates, bond DP, site DP and 
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Figure 5. Confirmation of universality. The figure shows the 3-bit probabilities Sc{b) 
with b ^ 2,4, 8, 16, 32, 64 for various well-known lattice models of DP, namely directed 
bond percolation, the contact process, directed site percolation and Wolframs W18 
stochastic cellular automaton. In addition, the figure shows data for the Domany- 
Kinzel model with optimized parameters (see text). As can be seen, all estimates tend 
to the same universal values as 6 — oo. 



Wolframs rule W18 (see [10] and references therein). As shown in Fig. |5l these models 
produce different values for Sc{b) but they converge towards the same value as the block 
size is taken to infinity, confirming the universality hypothesis. Surprisingly the contact 
process and bond percolation produce almost identical results, indicating that they are 
probably characterized by very similar non-universal short-range correlations. 

Moreover, we studied the Domany-Kinzel (DK) model [18] along its critical line 
parametrized by pi and p2, searching for a point where the deviation for small block 
sizes is minimal (see Fig. E]). Without going into detail we would like to report that this 
optimal point is located at pi = 0.58640(2) and p2 = 0.95923(2). We expect that this 
transition point of the critical line is characterized by minimal corrections to scaling. 

4. Parity-conserving universality class 

Parity conserving (PC) models belong to a universality class of absorbing state 
transitions different from DP. The defining feature of this class of systems is that the 
dynamics conserves the number of particles modulo 2 which is referred to as the parity 
of the system. There are several variants of parity conserving models which are all 
characterized by the critical exponents [TD] 



P = 0.92(2), = 1.83(3), 



3.22(6). 



(24) 
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Figure 6. Two-block probabilities for the parity-conserving class. Left panel: Sc{b,t) 
as a function of 1 /t for various block sizes b measured in the PC model defined in (^5]) . 
Right panel: Extrapolated values Sc{b) = limt_>oo Sc{b,t) plotted against 1/6. 



To apply the block renormalization procedure for the PC class, one needs to make 
sure that the renormalization scheme also respects the additional microscopic symmetry 
present in the system, i.e. the conservation of parity. This can be ensured if the 
reduction of a block of length 6 to a single site keeps the parity of that block unchanged. 
This is not possible with a two-state model, where a bit obtained by reducing a block 
does not carry any information about the parity. Instead one needs three states, namely, 
the inactive state and two active states indicating whether a block contains an even or 
odd number of particles. 

In this spirit we devise a parity-conserving three state model with an absorbing 
phase transition. As will be discussed below, this model renders itself easily to 
the scheme of block renormalization without hampering the parity of the original 
configuration. 

The model is defined on a one- dimensional periodic lattice, where each lattice site i 
is either empty or occupied by an odd or even number of particles, corresponding to the 
state variables Sj = 0, 1 and 2, respectively. The system evolves in time according to 
the following dynamical rules. 




with rate p, 

with rate 1 — p. (25) 

This dynamics either adds or removes a pair particles, keeping the parity of the 
system conserved. The only absorbing configuration is 000... since any occupied site 
is considered as active. The order parameter p is the density of active sites, not the 
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Figure 7. Extrapolated three-block probabilities Sc{b) plotted against 1/b. 
Extrapolating the lines to 6 — >■ 00 one obtains the numerical estimates of the universal 
numbers Sc which are listed in Table 1. 



density of particles in the system. This system shows an absorbing state transition 
belonging to the PC class defined by the critical exponents f l24p at the critical point 

= 0.4017(1) . (26) 

The parity-conserving variant of the block renormalization scheme can be applied to 
this model as follows. A block of length b is replaced by a single site with the state 

if the block is empty, 

1 if the total parity with the block is odd, 

2 if the total parity with the block is even. 

This choice makes sure that the parity of the renormalized configuration is same as that 
of the original one. Of course, as in the previous case, the absorbing configuration of 
the original lattice still remains absorbing under this renormalization scheme. 

Measurement of the conditional probabilities 

Even in the parity-conserving case one can directly apply the analysis of section [2] in 
order to determine the universal numbers Sc- However, as an active site i can be in two 
possible states Sj = 1, 2, we use the subscripts '0' and 'a' to denote inactive and active 
sites. For example, the renormalized density of active sites is -Pa(^) while S'aa denotes 
the conditional probability of finding two consecutive active blocks and so on. It is 
then straight-forward to prove that the exact expressions fl20p and (12 ip of the two-bit 
probabilities in terms of the critical exponents generalized to the parity-conserving case 
have the same form in the new notation, i.e. 

So. = Sao = 1 - 2-^^/^^ , 5aa = "1 + 2^-'^/^^ . (27) 
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Using the estimates fl2^ for the critical exponents of the PC class, we get 

^Oa = 5aO = 0.294(3) , 5aa = 0.411(2). (28) 

In order to verify this prediction, we simulated the PC model defined in (125|) with a 
system size L = 10^. Starting from a fully active configuration, the system evolves at the 
critical point ( 126|) for 8x10^ Monte Carlo cycles. In certain time intervals, b consecutive 
sites of the actual configuration of L sites are reduced to a single site by applying the 
renormalization scheme described above. This block-renormalized configuration yields a 
string of L/b sites. Again the conditional probabilities Sc{b, t) are measured by scanning 
this string and averaging over many runs for block sizes 6 = 1, 2, 4, 8, 16, 32. 

Following the same procedure as in the previous section, the left panel of Fig. E] 
shows the two-block quantities Sc{b, t) as functions of 1/t. The asymptotic values 5'aa(&) 
and S'ao(fe) obtained by extrapolating these curves to t — )■ oo are then plotted against 
l/b in the right panel of Fig. |6l Finally, the universal numbers Sc = limb_j.oo Sc{b) are 
estimated by another linear extrapolation of these curves. As can be seen, they are in 
excellent agreement with the predicted values (128|) . 

We have also measured the conditional probabilities Sc{b) for three adjacent blocks, 
for which no theoretical prediction is available. Fig. [7] shows plot of Sc{b) vs l/b for 
different patterns c, where we distinguished only inactive and active blocks. Note that 
in this case there are only five independent curves because of the spatial reflection 
symmetry Sa^Qo{b,t) = SoQa.{b,t) and S'aao(^5^) = SQ^s.{b,t). Clearly, these values approach 
a well defined limit Sc which is estimated by linear extrapolation as before. A summary 
of all estimates is given in Table 1. 

5. Summary 

In this paper we have introduced a simple renormalization scheme suitable for systems 
with absorbing states. For two-state systems the renormalization step involves a logical 
OR of the elementary degrees of freedom. The scheme can be applied in any space 
dimension to spatial as well as spatio-temporal configurations. 

The present work is restricted to the analysis of spatial configurations in 1+1 
dimensions. It is shown that in a critical process with homogeneous initial conditions 
all block probabilities Pc{b,t) for non- vanishing configurations decay asymptotically 
as differing only by their proportionality constants, which means that their ratios 
will tend to certain constants as t — ?■ oo. This allows us to define (non-universal) 
ratios Sc{b) = limt^oo Pc{b,t) / {1 — Po...oib,t)) which depend only on the block size. As 
the renormalization procedure has the tendency to suppress non-universal microscopic 
details, these quantities converge to universal values Sc in the limit of large block sizes 
6 — )■ oo. It is important to keep in mind that the two limits, t — oo and 6 — )■ oo, do not 
commute. 
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Table 1. Estimates of the universal conditional probabilities Sc for directed 
percolation and the parity-conserving universality class, c' denotes the spatially 
reflected bit pattern. The 2-bit values were obtained by using the scaling relations 
(Uni) I (EH) and (pS)) . the other values by numerical simulations. 

Using the examples of directed percolation and the parity-conserving class, we have 
demonstrated that the universal numbers are as robust as the critical exponent and 
therefore characterize the underlying universality class of the process like a fingerprint. 
The two-bit quantities turn out to depend on the critical exponents by exact scaling 
relations, providing a numerical test for the accuracy of the method. The three- and four- 
bit quantities, however, render larger sets of universal numbers which are independent 
of the critical exponents. By comparing the numerical estimates for a given model with 
the values listed in Table 1, this method provides a very sensitive test to which of the 
universality classes the phase transition belongs. As renormalization drives a system 
away from criticality, we note that it is crucial for the method to use a very accurate 
estimate of the critical point. 

The proposed method could shed new light onto nonequilibrium processes whose 
critical behavior is still debated in the literature. Examples include the pair contact 
process with diffusion and conserved sandpile models, which will be analyzed in a 
forthcoming paper. 
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